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Abstract 

We present an accurate computation of the ^-factors of the hyperhne states of the hydrogen 
molecular ion H^. The results are in good agreement with previous experiments, and can be 
tested further by rf spectroscopy. Their implication for high-precision two-photon vibrational 
spectroscopy of is also discussed. It is found that the most intense hyperhne components of 

j> ■ two-photon lines benefit from a very small Zeeman splitting. 
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I. INTRODUCTION 

In we have studied the spectrum of two-photon ro- vibrational transitions in the 
hydrogen molecular ion H^. The precise measurement (using a Doppler-free excitation 
geometry) of the frequency of such transitions in a rf trap constitutes a promising new 
method for determination of the electron-to-proton mass ratio m e /m p 2, Our estimate 
of transition rates with present-day experimental parameters has shown the feasibility of 
such an experiment. In order to assess the metrological merits of the method, it is essential 
to evaluate the order of magnitude of systematic effects; if one of them turns out to be a 
limiting factor, it is desirable to calculate it precisely in order to subtract its effects from 
the experimental data. 

Among the systematic effects, the ac Stark shifts due to the optical fields, blackbody 
radiation and rf trap potential are expected to be small, due to the small dynamic polariz- 
ability of the ion [4|. The second-order Doppler effect results in a shift and broadening 
of the line by about 10 kHz in a typical Paul trap, but can be reduced by sympathetic ion 
cooling [5] . The Zeeman shift is one of the most important effects remaining to be investi- 
gated, especially if circular polarization is used, resulting in the selection rule AMj = ±2 
for two-photon transitions. In the low magnetic field regime, the hyperfine structure has to 
be taken into account. The g-factors of the hyperfine levels have been calculated long ago [f| 
but with an accuracy limited by an imperfect knowledge of the hyperfine structure. A mea- 
surement of ^-factors ratios has been used by Richardson et al. to extract improved values 
of the hyperfine Hamiltonian coefficients [3] • We have recently determined these coefficients 
ab initio with an improved relative accuracy of 0(a 2 ), corresponding to the limit of the 
Breit-Pauli Hamiltonian 8]. The main aim of the present work is to calculate the ^-factors 
with much better accuracy than obtained so far. This high accuracy is, in fact, not needed 
for optical spectroscopy experiments, where the magnetic field is usually not controlled very 
precisely, so that the uncertainty on the Zeeman shift will be limited by the uncertainty 
on the magnetic field and its eventual fluctuations. However, we think it is worthwhile to 
present these precise values, because they can be tested in rf spectroscopy experiments of 
the same type as described in [7], which would provide a good test of hyperfine structure 
calculations. 
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II. ZEEMAN HAMILTON! AN 



Neglecting relativistic and radiative corrections, the linear part of the Hamiltonian de- 
scribing interaction of a ion with a magnetic field is given by 

H z = 9e^B S e • B - g p n P I • B + /i B L e • B - fi p (L x +L 2 ) • B, (1) 

where g e and g v are respectively the electron and proton g-factors, /ig = e/2m e is the 
Bohr magneton, and L e ,L l5 L 2 are the orbital momenta of the electron and both protons 
in the center-of-mass frame. S e and I are the electron and nuclear spins (see Ref. l|). The 
magnetic field is assumed to be oriented along Oz. In the ground electronic state lscx 9 , 
the main contribution to the ^-factor comes from the first term, i.e. the electron magnetic 
moment, while the other terms are about a factor of 1000 smaller. For the third term this 
can be understood by noting that the rotation velocities of the electron and protons are of 
the same order; hence the terms in L e and L1 + L2 are of the same order. 

In low magnetic fields, the hyperfine structure has to be taken into account, and the 
Zeeman Hamiltonian can be written using the Lande factor (or ^-factor) of the hyperfine 
level under consideration: 

H z = gj^B J B (2) 

where J is the total angular momentum; note that the scalar gj becomes a tensor if relativis- 
tic and radiative corrections are taken into account, in contradistinction with the atomic 
case, due to the lack of central symmetry of the potential 9|. 

Let us briefly recall the structure of hyperfine levels (for more details, see For even 
values of L, the total nuclear spin / is zero the total angular momentum is J = L + S e , and 
each ro-vibrational level (v , L) is split into two hyperfine levels (v,L,J) with J = L±l/2. 
For odd values of L, I is equal to one and the coupling scheme is as follows: the total spin 
F = S e + 1 (F = l/2 or 3/2) which is not an exact quantum number, and the total angular 
momentum J = L + F. J can take the values L±3/2, L±l/2. If J = L±3/2, the hyperfine 
state is a pure F = 3/2 state \v, L, S e = 1/2, 1= 1, F = 3/2, J = L±3/2,Mj). If J = L ± 1/2, 
the hyperfine state is a linear combination of F = l/2 and F = 3/2 states: 

I v, L, S e , I, F, J=LU/2, Mj) = Cf \ v, L, 1/2, 1, 1/2, LB/2, Mj)+Cf \ v, L, 1/2, 1, 3/2, LU/2, Mj). 

(3) 
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In the following, we evaluate the contribution of each term in equation (pQ) to the g-factor 
of the pure states \v, L, S e , I, F, J, Mj), noted as gj{v , L, F, J) (if L is even, these notations 
are replaced by \v, L, S e , J) and gj(v, L, J), respectively). The effect of state mixing will be 
addressed in Sec. IVHI The g-factor is divided into three contributions: 

gj(v, L, F, J) = 9l (L, F , J) + g 2 (L, F, J) + g 3 (y, L, F, J) (4) 

where g\ is the contribution from the first term in equation (pQ), gi is the contribution from 
the second term, and g% the contribution from the last two terms. Note that the first two 
quantities do not depend on the vibrational quantum number v, because the Hamiltonian 
only involves spin operators. In contradistinction, the last term contains orbital momentum 
operators acting on the orbital wave function, which introduces a slight dependence on v. 
The standard angular algebra procedures used below can be found in many textbooks, 



e.g. (l0| 



III. CONTRIBUTION OF THE ELECTRON SPIN 

Here, we evaluate the contribution to the ^-factor coming from the first term in equa- 
tion ([T]). The Zeeman shift of a given hyperfine level \v, L, S e , I, F, J, Mj) due to this term 
is given by 

AE ± = g ef i B B (v, L, S e7 J, F, J, Mj\S z e \v, L, S e , I, F, J, Mj) . (5) 

From now on the dependence in v will not be noted, since SI does not act on the orbital 
wave function. Application of the Wigner-Eckart theorem to the vector operator S e yields 

(L,S e ,I,F,J,Mj\S*\L,S e ,I,F,J,Mj) = -=L= (L, S e , I, F, J\\S e \\L, S e , I, F, J) (J1MjO\JMj) 

}Ij (L,S e ,I,F,J\\Se\\L,S e ,I,F,J). 



/ J(J+1)(2J+1) 

(6) 

The contribution to the ^-factor of the hyperfine level under study is then: 

u t? T \ (L,S e ,I,F,J\\S e \\L,S e ,I,F,J) 
g 1 (L,F,J)=g e " = (7) 

rj{J + l){2J + l) 



We now consider separately the cases of even and odd L. 
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A. Even values of L 



In this case the intermediate angular momentum F is irrelevant since I = 0, and one 
directly has J = L + S e . The reduced matrix element appearing in (j7|) is then obtained as 



(LS e J\\S e \\LS e J) = (S e \\S e \\S e ) (-1) 



J+L+S e + 1 



(27+1) 



S e 1 S e 
J L J 



(8) 



where (S , e ||S' e ||S' e ) 
values of 7: 



S e (S e +l)(2S e + l) 



9l (L,L+l/2) 
9l (L,L-l/2) 



3/2 . One gets g\(L, 7) for the two possible 

9c 



2L + V 

9e 

2L + 1 



(9) 



B. Odd values of L 



In this case the total spin of nuclei is one: 1 = 1. The reduced matrix element appearing 
in ([7]) is obtained in two steps: 



F I F 



(S e ,I,F\\S e \\S e ,I,F) = (S e \\S e \\S e ) (2F + 1) 

(L,S e ,7,F,7||S e ||L,S e ,7,F,7> = (S e , 7, F||S e ||S e , 7, F> ( 2 J+1) 
= J- (_i)J+z+i/2 (2F + 1)(27 + 1 



(10) 



F 1 F 
J L J 





1 S e 




F 1 F 






1! 




\' 


7 F 

> 







From this we get the factors g\ (L, F, 7) for all hyperfine levels: 

9e 1 



9l (L,l/2,L+l/2) : 
9l (L,l/2,L-l/2) 
5l (L,3/2,L+3/2) = /y , 
5l (L,3/2,L+l/2) : 
9l (L,3/2,L-l/2) 
9l (L,3/2,L-3/2) = -g e 



3 2L + 1 
ge 1 
3 2L + 1 
1 



2L + 3 
9e 2L + 9 



(12) 



3 (2L + l)(2L + 3) 
9e 2L-7 
' 3 (2L-1)(2L + 1) 
1 



2L - 1 
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IV. CONTRIBUTION OF THE NUCLEAR SPIN 



In the same way, we evaluate the contribution from the second term in equation ([T]), 
which is nonzero only for odd values of L. Similarly to equation (j7|), the ^-factor associated 
with this term is 

(T T? J\ m e ( L ^ S e,I,F,J\\I\\L,S e ,I,F,J) 

g 2 (L, F, J) = -g p _ (13) 

m P sj J{J +1)(2 J +1) 

The reduced matrix element appearing in the above equation is obtained similarly to equa- 
tion (HI]): 



(L, S e , I, F, J\\I\\L, S e , I, F, J) = y/6 (-1)^+V2 ( 2 F + 1)(2 J + 1) 



From this we deduce the factors g%{L, F, J) for all hyperfine levels: 

, , m P 4 1 

ff2(i ,l/2,L + l/2) = - 9p -i.-_ , 

771 A \ 

g 2 (L,l/2,L-l/2) = g P ^-- 



I 1 I F 1 F 
F S e F\\ J L J 



(14) 



m p 3 2L + 1 



£ 2 (L,3/2,L + 3/2) = * 



m p 2L + 3 ' 

m e 2 2L + 9 
g2 (L,3/2,L + l/2) = - giJ -.- (2L+i)(2L + 3) . 

, , 2 2L — 7 

92 (L,3/2,L-l/2) = ft^-3 (2L _ 1)(2L+1) . 

?7Z 2 

<? 2 (L,3/2,L-3/2) = s P " 



(15) 



m p 2L — 1 



V. CONTRIBUTION OF THE ORBITAL MOMENTA 



The contribution to the ^-factor coming from the third and fourth terms in equation ([I]) 



is 



(u, L, ^e, 7, F, J||L e ||u, L, S e , I, F, J) m e (v, L, S e , I, F, J\\Li\\v, L, S e , I, F, J) 
g 3 {v,L,F,J) = 2 



'J(J + 1)(2 J +1) m v yJJ(J + l)(2J + l) 

(16) 

where we have used that (L\) = (L 2 ) due to the symmetry of with respect to the 
exchange of nuclei. The reduced matrix elements appearing in the above expression are 
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expressed as a function of reduced matrix elements involving only the orbital wave function: 

(v,L,S e ,I,F, J\\Li\\v,L,S e ,I,F, J) = {-l) J+L+F+1 (2J + 1) J L 1 L \ (v, L\\L t \\v, L) . 

\JFJ\ 

(17) 

This expression is valid both for even and odd values of L (in the first case F — S e — 1/2). 
We finally obtain the factors gs{v,L,F, J) for all hyperfine levels, expressed as a function 
of the orbital reduced matrix elements. For F = 1/2 levels (both with even and odd L) we 
have: 



ai (,,4l/U+l/^ ^ + 1) (| M ) 



g 3 (v,L,l/2,L-l/2) = JfF+J (\\L M \\) 

y L (zL + 1 ) 

and for F = 3/2 levels (appearing only if L is odd): 

*(^3/2,L + 3/2) = ^TT§^) (IM) 
53 (,,L,3/2,L+l/2) = . /r 2(2i ; 2 / +3L ~ 3) :(W) 



(18) 



2(2L 2 + L-4) 
/L(L + 1) 2L-1 (2L + 1 



where 



(19) 



ft(«,£.3/U-3/2) = ( ll^tll ) 

\ II tot 1 1 / ^/2L + l m p y/2L + l ' { ' 

The orbital matrix elements (v, L\\L e \\v, L) and (v, L\\Li\\v, L) have been calculated using 
the variational approach presented in |l|], f° r < t>, L < 4. We used basis lengths N = 
2000-3000, allowing a determination of the orbital momentum contribution g 3 with 10 -4 
relative accuracy. The results are shown in Table [I] 



VI. ROTATIONAL (/-FACTORS IN THE STRONG-FIELD REGIME 

These results allow a precise calculation of the rotational g-factor, which can be measured 
independently in a strong magnetic field. In this regime, the different angular momentum 
and spin vectors are decoupled, and the Zeeman Hamiltonian is usually written in the form 

H z = g e ^B S e B - g p fi p IB- g rot \i v L ■ B (21) 

7 



L 


V 


(u, L L e \\v, L) 


(v, L\\L\ \\v, L) 


(\\Ltot\\) 


9rot 


. /or i i 

V 1 


./or i 1 

V 1 























0.615e-04 


0.70708 


-0.7087e-03 


0.9201 




1 


0.686e-04 


0.70707 


-0.7015e-03 


0.9108 


1 


2 


0.763e-04 


0.70707 


-0.6938e-03 


0.9008 




3 


0.847e-04 


0.70706 


-0.6855e-03 


0.8900 




4 


0.937e-04 


0.70706 


-0.6764e-03 


0.8782 







1.069e-04 


1.22469 


-1.2271e-03 


0.9198 




1 


1.193e-04 


1.22469 


-1.2146e-03 


0.9105 


2 


2 


1.328e-04 


1.22468 


-1.2012e-03 


0.9004 




3 


1.473e-04 


1.22467 


-1.1867e-03 


0.8896 




4 


1.630e-04 


1.22466 


-1.1710e-03 


0.8778 







1.521e-04 


1.73197 


-1.7344e-03 


0.9193 




1 


1.698e-04 


1.73197 


-1.7167e-03 


0.9100 


3 


2 


1.889e-04 


1.73196 


-1.6977e-03 


0.8998 




3 


2.095e-04 


1.73195 


-1.6776e-03 


0.8889 




4 


2.318e-04 


1.73193 


-1.6547e-03 


0.8771 







1.980e-04 


2.23597 


-2.2375e-03 


0.9187 




1 


2.209e-04 


2.23596 


-2.2146e-03 


0.9092 


4 


2 


2.457e-04 


2.23595 


-2.1898e-03 


0.8991 




3 


2.725e-04 


2.23593 


-2.1629e-03 


0.8881 




4 


3.015e-04 


2.23592 


-2.1339e-03 


0.8761 



TABLE I: Reduced matrix elements of L e and Li (divided by y / 2L+T) for all ro-vibrational levels 
(v,L) of H2" with < v,L < 4. The deduced value of ( ||Ltoi|| ) (defined in equation ([20]) ) is 
indicated in the third column. The last column contains the values of the rotational g-factors 
deduced from equation (f22|) . 
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Comparing this expression with equation (pQ), one easily obtains the expression of the rota- 
tional (/-factor for a ro- vibrational level (v, L): 

9rot(v,L)= - 1 (2(v,L\\L l \\v,L)-^ (v, L\\L e \\v, L)) = ( , l|Ltot|l) 

^L(L + 1)(2L + 1) V m e / m e ^ L(L + I) 

(22) 

The values of rotational (/-factors are given in the last column of Table [B They are improved 
with respect to the previous calculation performed within the adiabatic approximation by 



Rebane and Zotev [131 ] . Note that the slight dependence of g rot on L is neglected in the 



approach followed in that paper 



Loch et al. 11[ measured the rotational (/-factor of , averaged over the vibrational 



levels v = 4 — 6 and the rotational levels L — 1 — 3, to g rot = 0.920(40). Extending our 



calcu 



ations to the v = 5 — 6 levels, and assuming the vibrational state populations reported 



in 11]. wc find g rot = 0.8688. in disagreement with the experimental value by 1.28 o. 



VII. EFFECT OF THE STATE MIXING 

We now have to take into account the fact that for odd L, some of the hyperfine eigenstates 
are not pure states with a well-defined value of F, but linear combinations of F — 1/2 and 
F = 3/2 states (see equation (j3J)). The various contributions to the (/-factor (denoted gj 
from now on) are changed in the following way: 

gj(v, L, F , J) = ~ 9l (v, L, F , J) + g 2 (v, L, F, J) + g 3 (v, L, F, J), (23) 

with 

Uv,L,F,L±l/2) = (C 1 ± ) 2 (/ 1 (L,l/2,L±l/2) + (C 3 ± ) 2 (/ 1 (L,3/2,L±l/2) 

^ . _ ± _ ± (L, S e , 1, 1/ 2, L±l/2||S e ||L, S e , 1, 3/2, L±l/2) 

T^i </e 



'J(J+1)(2J + 1) 

g 2 (v,L,F,L±l/2) = (C , 1 ± ) 2 (/ 2 (L,l/2,L±l/2) + (C 3 ± ) 2 (/ 2 (L,3/2,L±l/2) 

„ ± „ ± m e (L, g e , I 1/2, L±1/2||/||L, ff e , J, 3/2, L±l/2) 
m v Jj(J + 1)(2J+ 1) 



L, F, L±l/2) = (Cf) 2 (/ 3 (^, 1/2, L±l/2) + (Ct) 2 g 3 (v, L, 3/2, L±l/2). 

(24) 

The contribution g 3 does not contain any crossed terms, because the corresponding Hamil- 
tonian acts only on orbital variables, and cannot couple F — 1/2 states to F = 3/2 states. 



The mixing coefficients Cf, Cf are given in [12J. It only remains to calculate the crossed 



reduced matrix elements appearing in the above expression. We have 

{L,S e ,I,F,J\\S e \\L,S e ,I,F',J) = a (_i)W2^+3/2 ^ 1 Se ) {Se \\S e \\S e 

\F' I F 

(L, S e , I, F, J\\ / \\L, S e , I, F', J) = a (_l)'+*+s/» ) 1 1 1 / 



(25) 



with 



F' S P . F 



F 1 F' 



a = \l(2F + l)(2F' + l)(2J+l){ }. (26) 

J L J 

From these expressions one may finally obtain 

Js Jul + 1) 

(L, 1/2,1, 1/2, L±l/2\\S e \\L, 1/2,1, 3/2, L±l/2) = \- V _— 

o v^-^ + J- ('27) 
,/L(L + 1) 

(L, 1/2,1, 1/2, L±1/2\\I\\L, 1/2,1, 3/2, L±l/2> = ^=^ - 

VIII. FINAL RESULTS AND DISCUSSION 

Our results are summarized in Tables [TTl IHH and IIVI where all contributions are summed 
up in order to obtain the (/-factors for all hyperfine levels. The accuracy of this calculation 
can be limited by several factors. 

The first one is the variational calculation of the orbital momentum reduced matrix 
elements appearing in expression ( JT6l) of g%, the relative accuracy of which is 10 -4 . This 
results in a numerical uncertainty of less than 10~ 7 on the final (/-factor values. 

The second limitation comes from the coefficients of the hyperfine Hamiltonian, which 
are known with a relative accuracy 0(a 2 ) ~5 x 10~ 5 . The mixing coefficients defined in (J3]) 
are affected by this uncertainty, resulting in an uncertainty comprised between 5xl0 -7 and 
5xl0 -6 for the final values of g(L, F, J), depending on the magnitude of the corrections due 
to state mixing. This uncertainty does not affect the hyperfine states which are pure states 
(i.e. all states of even L, and states of odd L with F = 3/2 and J = L±3/2). 

Finally, and most importantly, relativistic and radiative corrections, as considered by 



Hegstrom in the strong- field regime in [9], are not included. This limits the relative accuracy 
to 0(a 2 )~5xl(T 5 . 
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L 


V 


J 


9i{L,J)/g e 


93{v,L,J) 


gj(v,L,J) 







1/2 


+1 





2.0023193 







3/2 


-1/5 


-6.011e-04 


-0.4010650 






5/2 


+1/5 


-4.008e-04 


0.4000631 




1 


3/2 


-1/5 


-5.950e-04 


-0.4010589 






5/2 


+1/5 


-3.967e-04 


0.4000672 


2 


2 


3/2 


-1/5 


-5.885e-04 


-0.4010523 






5/2 


+1/5 


-3.923e-04 


0.4000716 




3 


3/2 


-1/5 


-5.814e-04 


-0.4010452 






5/2 


+1/5 


-3.876e-04 


0.4000763 




4 


3/2 


-1/5 


-5.737e-04 


-0.4010375 






5/2 


+1/5 


-3.824e-04 


0.4000814 







7/2 


-1/9 


-5.559e-04 


-0.2230358 






9/2 


+1/9 


_4.447e-04 


0.2220352 




1 


7/2 


-1/9 


-5.502e-04 


-0.2230301 






9/2 


+1/9 


-4.402e-04 


0.2220398 


4 


2 


7/2 


-1/9 


-5.441e-04 


-0.2230240 






9/2 


+1/9 


-4.352e-04 


0.2220447 




3 


7/2 


-1/9 


-5.374e-04 


-0.2230173 






9/2 


+1/9 


-4.299e-04 


0.2220500 




4 


7/2 


-1/9 


-5.302e-04 


-0.2230101 






9/2 


+1/9 


-4.241e-04 


0.2220558 



TABLE II: g-factors of all hyperfine states for the ro-vibrational levels (v, L) with L = 0, 2, 4 and 
< v < 4. The fourth column is the contribution from the electron spin magnetic moment, the fifth 
one is the contribution from orbital momenta, and the last one is the total value of gj(v, L, J). All 
digits are converged. The relative theoretical accuracy on all contributions (and on gj) is 0(a 2 ). 

To sum it up, the accuracy of our results is 0(a 2 ) ~5x 10~ 5 . Inclusion of the leading- 
order QED corrections would improve it to 0(a 3 ) ~3x 1CT 7 . Such accuracy would appear 
to be hard to meet in rf experiments with weak magnetic fields, (/-factor measurements 
on with an accuracy in the ppm range were achieved in strong magnetic fields, using 
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V 


F 


J 


9l/9e 


m p 92 
mC ^ 


93 


gj(v,L,F, J) 


91 he 


mp 92 
e 


93 


gj(v,L,F, J) 




1/2 


1/2 


+ 1/9 


+4/9 


-6.682c- 


04 


+0.2231638 


+0.0629313 


+0.3932396 


-6.667c- 


04 


+0.1265381 






3/2 


-1/9 


—4/9 


-3.341c- 


04 


—0.2241660 


—0.1172298 


—0.4506607 


-3.340c- 


04 


— 0.2364365 





3/2 


1/2 


+5/9 


— 10/9 


+3. 341c- 


04 


+ 1.1093536 


+0.6037354 


— 1.0599063 


+3.326e- 


04 


+ 1.2059793 






3/2 


+11/45 


— 22/45 


-1.336c- 


04 


+0.4878350 


+0.2505632 


—0.4826727 


-1.337c- 


04 


+0.5001054 






of z 


+ l/o 


— ^/ o 


-2.005c- 


04 


+U. oyyU4DO 














1/2 


1/2 


+ 1/9 


+4/9 


-6.614e- 


04 


+0.2231705 


+0.0648181 


+0.3953621 


-6.600e- 


04 


+0.1303292 






3/2 


-1/9 


—4/9 


-3.307e- 


04 


—0.2241627 


—0. 1170219 


—0.4504461 


-3.307c- 


04 


— 0.2360162 


1 


3/2 


1/2 


+5/9 


— 10/9 


+3.307e- 


04 


+1.1093503 


+0.6018486 


— 1.0620288 


+3.293e- 


04 


+ 1.2021916 






3/2 


+11/45 


— 22/45 


-1.323c- 


04 


+0.4878363 


+0.2503552 


—0.4828872 


-1.323c- 


04 


+0.4996898 






5/2 


+ 1/5 


— 2/5 


-1.984e- 


04 


+0. 39904OO 










1 A OAAAviO/^ 

+ U.o99U4od 




1/2 


1/2 


+ 1/9 


+4/9 


-6.541e- 


04 


+0.2231778 


+0.0667416 


+0.3975158 


-6.529e- 


04 


+ 0.1341944 






3/2 


-1/9 


—4/9 


-3.271c- 


04 


—0.2241590 


—0.1168072 


—0.4502248 


-3.270c- 


04 


— 0.2355819 


2 


3/2 


1/2 


+5/9 


— 10/9 


+3. 271c- 


04 


+1.1093466 


+0.5999250 


— 1.0641825 


+3.258e- 


04 


+ 1.1983300 






3/2 


+11/45 


— 22/45 


-1.308c- 


04 


+0.4878378 


+0.2501405 


—0.4831085 


-1.309c- 


04 


+0.4992607 






5/2 


+ 1/5 


— 2/5 


-1.962e- 


04 


1 A OAAACAO 

+0.3990b0o 










1 A OAAACAO 




1/2 


1/2 


+ 1/9 


+4/9 


-6.463e- 


04 


+0.2231857 


+0.0687069 


+0.3997057 


-6.451e- 


04 


+0.1381439 






3/2 


-1/9 


—4/9 


-3.231c- 


04 


—0.2241551 


—0.1165853 


—0.4499964 


-3.231e- 


04 


— 0.2351330 


3 


3/2 


1/2 


+5/9 


— 10/9 


+3. 231c- 


04 


+1.1093427 


+0.5979598 


— 1.0663724 


+3.220e- 


04 


+ 1.1943845 






3/2 


+ 11/45 


— 22/45 


-1.293c- 


04 


+0.4878393 


+0.2499186 


—0.4833369 


-1.293c- 


04 


+0.4988172 






5/2 


+ 1/5 


-2/5 


-1.939c- 


04 


+0.3990532 










+0.3990532 




1/2 


1/2 


+ 1/9 


+4/9 


-6.377c- 


04 


+0.2231942 


+0.0707213 


+0.4019394 


-6.367e- 


04 


+0.1421927 






3/2 


-1/9 


-4/9 


-3.189c- 


04 


-0.2241508 


-0.1163546 


-0.4497592 


-3.188c- 


04 


-0.2346661 


4 


3/2 


1/2 


+5/9 


-10/9 


+3. 189c- 


04 


+1.1093384 


+0.5959453 


-1.0686061 


+3.179e- 


04 


+ 1.1903399 






3/2 


+ 11/45 


-22/45 


-1.275c- 


04 


+0.4878410 


+0.2496879 


-0.4835741 


-1.276c- 


04 


+0.4983563 






5/2 


+ 1/5 


-2/5 


-1.913e- 


04 


+0.3990557 










+0.3990557 



TABLE III: ^-factors of all hyperfine states for the ro-vibrational levels (v, V) with L = 1 and 
< v < 4. The fourth column is the contribution from the electron spin magnetic moment, the 
fifth one is the contribution from the nuclear spin magnetic moment, and the sixth one is the 
contribution from orbital momenta. All three terms are calculated without taking state mixing 
into account. The seventh column is the sum of these contributions. Columns 8 to 10 are the 
corrected values of the three contributions obtained by taking state mixing into account. The last 
column is the final value of the g-factor. All digits are converged. The relative theoretical accuracy 
for all contributions is 0(a 2 ). 



spin-dependent charge-exchange techniques 11]. 

We have nevertheless given the (7-factors values with 7 digits. Although all digits are not 
significant, this is helpful for understanding the order of magnitude of various effects, such as 
the variation as a function of v and L, or the importance of corrections due to state mixing. 
For example, it can be seen that state mixing acts on the value of the orbital contribution (73 
at the level of 10~ 6 at most, so this effect may be neglected at the present level of accuracy. 

To our knowledge, the magnetic moments of the H^" hyperfine states in a weak magnetic 
field have been investigated only by Richardson, Jefferts and Dehmelt in 1968 [7j]. They 
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V 


F 


J 


ail an 


m p j?2 
e 9p 


93 


3j(v,L,F, J) 


ffl/9e 


m P 92 
m e g v 


33 


gj(v,L,F,J) 




1/2 


5/2 


+1/21 


+4/21 


-5.722e 


-04 


+0.0953558 


+0.0264523 


+0.1686545 


-5.719c 


-04 


+0.0529071 






7/2 


— 1/21 


—4/21 


-4.292o 


-04 


—0.0963571 


—0.0566485 


— 0.1996746 


-4.291e 


-04 


—0.1144649 





3/2 


3/2 


-1/5 


+2/5 


-8.011c 


-04 


—0.4000481 










— 0.4000481 






5/2 


+ 1/105 


— 2/105 


-4.864e 


-04 


+0.0185254 


+0.0306905 


+0.0027741 


-4.867e 


-04 


+0.0609740 






7/2 


+5/63 


— 10/63 


-3.815e 


-04 


+0.1580499 


+0.0883945 


— 0.1495318 


-3.816c 


-04 


+0.1761576 






9/2 


+ 1/9 


— 2/9 


-3.338c 


-04 


+0.2214701 










+ 0.2214701 




1/2 


5/2 


+1/21 


+4/21 


-5.664e 


-04 


+0.0953616 


+0.0272667 


+0.1695181 


-5.661c 


-04 


+0.0545462 






7/2 


— 1/21 


—4/21 


-4.248c 


-04 


—0.0963528 


—0.0563438 


— 0.1993583 


-4.247c 


-04 


—0.1138494 


1 


3/2 


3/2 


-1/5 


+2/5 


-7.929c 


-04 


—0.4000400 










— 0.4000403 






5/2 


+ 1/105 


— 2/105 


-4.814e 


-04 


+0.0185303 


+0.0298762 


+0.0019104 


-4.817c 


-04 


+0.0593457 






7/2 


+5/63 


— 10/63 


-3.776c 


-04 


1 A 1 CO A COO 

+U. looUooo 


i a ncununc 


A 1 A AO A O A 

— U. 1498480 


-3.777c 


-04 


+ U. 1 ( OOOU4 






9/2 


+ 1/9 


—2/9 


-3.304e 


-04 


+0.2214735 










+ 0.2214735 




1/2 


5/2 


+1/21 


+4/21 


-5.601e 


-04 


i a nncosw 
+U.UyDo6 / y 


1 A AAOAAAA 


1 A 1 *7A'JAA'7 

+o.i i Ooyy i 


-5.598e 


-04 


+ U.(JOOZZo4 






7/2 


— 1/21 


—4/21 


-4.201e 


-04 


—0.0963480 


—0.0560290 


— 0.1990321 


-4.200c 


-04 


—0.1132134 


2 


3/2 


3/2 


-1/5 


+2/5 


-7.841c 


-04 


—0.4000311 










— 0.4000313 






5/2 


+ 1/105 


— 2/105 


-4.761c 


-04 


+0.0185357 


+0.0290430 


+0.0010288 


-4.763c 


-04 


+0.0576802 






7/2 


+5/63 


— 10/63 


-3.734c 


-04 


+0.1580580 


+0.0877750 


—0.1501742 


-3.735c 


-04 


+0.1749233 






9/2 


+ 1/9 


— 2/9 


-3.267c 


-04 


+0.2214772 










+0.2214772 




1/2 


5/2 


+1/21 


+4/21 


-5.533e 


-04 


+0.0953747 


+0.0289535 


+0.1713009 


-5.530c 


-04 


+0.0579423 






7/2 


— 1/21 


—4/21 


-4.150e 


-04 


—0.0963429 


—0.0557027 


— 0.1986944 


-4.149c 


-04 


— 0.1125539 


3 


3/2 


3/2 


-1/5 


+2/5 


-7.746c 


-04 


—0.4000216 










— 0.4000217 






5/2 


+ 1/105 


— 2/105 


-4.703e 


-04 


+0.0185415 


+0.0281893 


+0.0001276 


-4.705c 


-04 


+0.0559739 






7/2 


+5/63 


— 10/63 


-3.689c 


-04 


+0.1580625 


+0.0874487 


— 0. 1505119 


-3.689c 


-04 


+0.1742735 






9/2 


+ 1/9 


-2/9 


-3.227c 


-04 


+0. 2214812 










i n 001 i 
+ U.^zl4»lz 




1/2 


5/2 


+ 1/21 


+4/21 


-5.459e 


-04 


+0.0953821 


+0.0298308 


+0.1722248 


-5.457c 


-04 


+0.0597089 






7/2 


-1/21 


-4/21 


-4.094c 


-04 


-0.0963374 


-0.0553637 


-0.1983442 


-4.094c 


-04 


-0.1118686 


4 


3/2 


3/2 


-1/5 


+2/5 


-7.643c 


-04 


-0.4000113 










-0.4000115 






5/2 


+ 1/105 


-2/105 


-4.640c 


-04 


+0.0185477 


+0.0273121 


-0.0007963 


-4.642c 


-04 


+0.0542209 






7/2 


+5/63 


-10/63 


-3.639c 


-04 


+0.1580674 


+0.0871097 


-0.1508622 


-3.640c 


-04 


+0.1735986 






9/2 


+ 1/9 


-2/9 


-3.184e 


-04 


+0.2214855 










+0.2214855 



TABLE IV: Same as Table 1111} for the rotational level L = 3. 



give a few ratios of g-factors between different hyperfine states, which we have reported in 
Table |V] together with the result of our calculation. Note that the experimental values are 
averaged over the vibrational states v > 4. In our evaluation, we have taken the vibrational 
states v = 5-8 into account, and assumed that their relative populations (determined by the 
creation process by electron impact ionization of H 2 at room temperature) are the same as 



measured in Ref. 



14| . Good agreement is obtained in all cases. 



Ratio 


Calculated 


Measured [7J 


gj{L = l,F = l/2,J=3/2)/gj(L=l,F=3/2,J=5/2) 
gj(L = l,F = 3/2,J = 3/2) / gj (L = 1, F = 3/2, J = 5/2) 
g j(L = 1, F = 1/2, J = 3/2) / gj (L = 3, F = 3/2, J = 9/2) 


0.5855 
1.2463 
1.0549 


0.584(3) 
1.241(6) 
1.051(5) 



TABLE V: ^-factor ratios of some hyperfine states of Hg" 
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We now use these results to evaluate the Zeeman shift and splitting of the two-photon 
transitions (v — 0, L) — > (v' — l,L) envisaged for high-precision spectroscopy of the Kj" ion. 
For illustration, we choose a magnetic field of the order of the earth field, B = 5x 10~ 5 T. 
As explained in the introduction, the accuracy of our calculation is more than sufficient 
for this purpose. The leading relativistic corrections which we have neglected correspond 
to a shift of order a 2 fj, B B ~ 35 Hz, well below the present goal accuracy of spectroscopy 
experiments. The frequency shift Av of a two-photon transition (v, L, F, J) — > (v', L, F', J') 
in the magnetic field B is 



2hAv 



M'j~gj{v'i L, F', J') - Mj~gj(v, L, F, J)} ^ B B (28) 



If circular polarization is used, the selection rule is M'j — Mj = 2. Mj can take the possible 
values — J, — J + 1...J — 2 so that the shift of the line centre (corresponding to Mj —— 1) is 



2hAv 



gj(v, L, F, J) + gj(v', L, F\ J') ^ B B (29) 



In the case of linear polarization, the selection rule is Mj — Mj = 0, and Mj can take the 
possible values — J, — J + 1...J. There is no global shift, and the Zeeman splitting between 
extreme values of Mj is 

hAv = [gj(v\ L, F', J') - gj(v, L, F, J)] Jfi B B (30) 

We have shown in []J that the most intense hyperfine components are those between pairs 
of homologous spin states, (F, J) — > (F,J), and only these components are considered in 
the following. In this case, we benefit from an almost complete cancellation (to 1 percent or 
better) between the ^-factors of the initial and final states, so that the Zeeman splitting is 
very small (compared to the global shift observed in circular polarization). Note that such 
cancellation will also take place for relativistic corrections, so that the theoretical uncertainty 
is also reduced. In estimating the uncertainties, we have assumed cancellation to 1 percent. 



The cases L = 1,2,3 are compared in Table IVTl for the cases of circular and linear 
excitation polarizations. When circular polarization is used, the two-photon transition lines 
are typically shifted by a few hundreds of kHz. This does not represent in itself a limitation 
of experimental accuracy, since it is possible take the average of measurements in a + and cr_ 
polarizations. However, it also means that magnetic field fluctuations of the order of 10 mG 
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L 


F 


J 


shift 


splitting 








(<T + polar.) 


(ir polar.) 





1/2 


1/2 


* 







1/2 


1/2 


* 


1327 






3/2 


-165 314(9) 


441 


1 


3/2 


1/2 


* 


1325 






3/2 


+349 834(19) 


436 






5/2 


+279 258(15) 


3.5(3) 


2 


1/2 


3/2 


-280 668(15) 


6.4(2) 




1/2 


5/2 


+279 971(15) 


7.2(3) 




1/2 


5/2 


+ 37 599(2) 


2868 






7/2 


- 79 889(4) 


1508 


3 


3/2 


3/2 


-279 956(15) 


8.2(2) 






5/2 


+ 42 101(2) 


2849 






7/2 


+123 065(7) 


1487 






9/2 


+154 989(8) 


10.7(3) 



TABLE VI: Zeeman shift and splitting of two-photon transition lines (v = 0, L, F, J) — > (v' = 
1,L,F,J) in a magnetic fields of 5 x 10 -5 T, in Hz. The fourth column is the predicted shift 
of the line center in the case of o+ excitation polarizations (equation (|29p ). The stars indicate 
transitions which are forbidden in circular polarization. An estimate of the theoretical uncertainty 
(corresponding to a relative accuracy of 0(a 2 )) is given. The last column is the Zeeman splitting, 
i.e. the frequency difference between extreme values of Mj, evaluated in the linear polarization case 
(equation (|30j) ). An estimate of the theoretical uncertainty (see text) is given when it is significant. 



result in a line broadening of order 1-10 kHz. If one wishes to improve the resolution beyond 
this limit, active control and stabilization of the magnetic field is required Together 
with larger transition probabilities as discussed in lj, this brings a strong argument in 
favor of using linear polarizations. Unfortunately, optical isolation of the laser source from 
feedback by the enhancement cavity is at present only possible with a polariser followed by 
a quarter- wave plate, which imposes working with circular polarizations \L\. One solution is 
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to add a transverse magnetic field, which must be sufficiently strong to separate the three 
components AMj = 0, ±2. It can be seen from Table [VTl that a field in the 10~ 5 -10~ 4 T 
range (depending on the transition) is enough to obtain a separation of the order of 100 
kHz, i.e. clearly resolved components assuming a linewidth of a few kHz [l]. 

Finally, the Zeeman splitting is extremely small (a few Hz) when the states involved in 
the two-photon transition are pure states, while it is of a few kHz in other cases. This makes 
such transitions especially attractive from the metrological point of view. 



IX. CONCLUSION 



We have obtained improved values of (/-factors of the hyperfine states of the hydrogen 
molecular ion, which are in good agreement with experiment. The achieved accuracy is 
0(a 2 ) ~ 5x 10~ 5 . The accuracy on the rotational (/-factors has also been improved by use 
of a variational method allowing to take the full three body dynamics into account. We 
have used these results to evaluate the Zeeman shift and splitting of several two-photon 
vibrational transition lines, and shown that transitions involving pure hyperfine states (i.e. 
all states of even L, and states of odd L with F = 3/2 and J = L±3/2) benefit from a very 
small Zeeman splitting. 
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